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Abstract 

The cross section for the removal of high-momentum protons from ^^O is cal- 
culated for high missing energies. The admixture of high-momentum nucleons 
in the ^^O ground state is obtained by calculating the single-hole spectral func- 
tion directly in the finite nucleus with the inclusion of short-range and tensor 
correlations induced by a realistic meson-exchange interaction. The presence 
of high-momentum nucleons in the transition to final states in ^^N at 60- 
100 MeV missing energy is converted to the coincidence cross section for the 
(e, e'p) reaction by including the coupling to the electromagnetic probe and 
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the final state interactions of the outgoing proton in the same way as in the 
standard analysis of the experimental data. Detectable cross sections for the 
removal of a single proton at these high missing energies are obtained which 
are considerably larger at higher missing momentum than the corresponding 
cross sections for the p-wave quasihole transitions. Cross sections for these 
quasihole transitions are compared with the most recent experimental data 
available. 

PACS Number(s): 21.10Jh.Jx, 21.30.+y, 24.10.Cn, 27.20.+n 
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I. INTRODUCTION 



Experimental progress in the exclusive (e, e'p) reaction in recent years has provided a 
clear picture of the limitations of the simple shell-model description of closed-shell nuclei. 
Of particular interest is the reduction of the single-particle (sp) strength for the removal of 
particles with valence hole quantum numbers with respect to the simple shell-model estimate 
which corresponds to a spectroscopic factor of 1 for such states. Typical experimental results 

for closed-shell nuclei exhibit reductions of about 30% to 45% for these spectroscopic 
factors. In the case of ^°^Pb, one obtains a spectroscopic factor for the transition to the 
ground state of '^^'^Tl of about 0.65 which is associated with the removal of a 3s| proton. An 
analysis which uses information obtained from elastic electron scattering, indicates that the 
total occupation number for this state is about 10% higher 0], corresponding to 0.75. This 
additional background strength should be present at higher missing energy and is presumed 
to be highly fragmented. The depletion of more deeply bound orbitals is expected to be 
somewhat less as suggested by theoretical considerations which also indicate that the 
strength in the background, outside the main peak, corresponds to about 10% (see also 0]). 

Recent experimental results for ^^O yield a combined quasihole strength for the p| 
and p| states corresponding to about 65% with the strength concentrated in one peak 
and the p| strength fragmented already over several peaks. Recent theoretical results yield 
about 76% for these p states without reproducing the fragmentation of the p| strength. 
This calculation includes the influence of both long-range correlations, associated with a 
large shell-model space, as well as short-range correlations. Although the inclusion of long- 
range correlations yields a good representation of the I = 2 strength, it fails to account for 
the presence of positive parity fragments below the first p| fragment. This suggests that 
additional improvement of the treatment of long-range correlations is indicated possibly 
including a correct treatment of the center-of-mass motion [0. The contribution to the 
depletion of the sp strength due to short-range correlations is typically about 10%. This 
result is obtained both in nuclear matter calculations, as reviewed in 0, and in calculations 
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directly for finite (medium-) heavy nuclei []7HT0|,|6|. Although the influence of long-range 
correlations on the distribution of the sp strength is substantial, it is clear that a sizable 
fraction of the missing sp strength is due to short-range effects. The experimental data [jl],|^ 
indicate that only about 70% of the expected protons in the nucleus has been detected in 
the energy and momentum domain studied so far. It is therefore important to establish 
precisely where the protons which have been admixed into the nuclear ground state due to 
short-range and tensor correlations, can be detected in the (e, e'p) reaction and with what 
cross section. 

The influence of short-range correlations on the presence of high-momentum components 
in finite (medium-) heavy nuclei has been calculated in [p|-[TO[]. In this work the spectral 



function for ^^O has been calculated from a realistic interaction without recourse to some 



form of local density approximation |irT| , p!2| . No substantial high-momentum components are 
obtained in |[^-[T0[| at small missing energy. With increasing missing energy, however, one 
recovers the high-momentum components which have been admixed into the ground state. 
The physics of these features can be traced back to the realization that the admixture of high- 
momenta requires the coupling to two-hole-one-particle (2hlp) states in the self-energy for 
a nucleon with high momentum. In nuclear matter the conservation of momentum requires 
the equality of the 2hlp momentum in the self-energy and the external high momentum. 
Since the two-hole state has a relatively small total pair momentum, one automatically needs 
an essentially equally large and opposite momentum for the intermediate one-particle state 
to fulfill momentum conservation. As a result, the relevant intermediate 2hlp states will lie 
at increasing excitation energy with increasing momentum. Considerations of this type are 



well known for nuclear matter (see e.g. '^^), but are approximately valid in finite nuclei as 



well. Recent experiments on ^osp]^ j-g g^j^^ i6q essentially confirm that the presence of 
high-momentum components in the quasihole states accounts for only a tiny fraction of the 
sp strength. 

The theoretical prediction concerning the presence of high-momentum components at 
high missing energy remains to be verified experimentally, however. In order to facilitate 



and support these efforts, the present work aims to combine the calculation of the spectral 
function at these energies with the description of both the electromagnetic vertex and fi- 
nal state interactions (FSI) in order to produce realistic estimates of the exclusive (e, e'p) 
cross section under experimental conditions possible at NIKHEF and Mainz. The impulse 
approximation has been adopted for the electromagnetic current operator, which describes 
the nonrelativistic reduction (up to fourth order in the inverse nucleon mass [|16|) of the 
coupling between the external virtual photon and single nucleons only. The treatment of 
FSI has been developed by the Pavia group |1T7|-^ (see also Ref. [^) and takes into ac- 
count the average complex optical potential the nucleon experiences on its way out of the 
nucleus. Other contributions to the exclusive (e, e'p) reaction are present in principle, such 
as two-step mechanisms in the final state or the decay of initial collective excitations in the 
target nucleus. However, by transferring sufficiently high energy uj to the target nucleus and 
by selecting typical kinematical conditions corresponding to the socalled quasielastic peak 
with uj = /2m {q the momentum transfer and m the nucleon mass), these contributions 
are suppressed. In these conditions, adopted in the most recent experiments, the direct 
knockout mechanism has been shown to be the dominant contribution and essentially 



corresponds to calculating the combined probability for exciting a correlated particle (which 
is ultimately detected) and a correlated hole such that energy and momentum are conserved 
but no further interaction of the particle with the hole is included. 

The calculation of the spectral function for ^^O is reviewed in Sec. II. Special attention 
is given to a separable representation of the spectral function which facilitates the practical 
implementation of the inclusion of FSI. In Sec. Ill the general formalism of the Distorted 
Wave Impulse Approximation (DWIA) is briefly reviewed. The influence of the FSI is studied 
in Sec. IV for the quasihole transitions for which data are available P,p!5[]. Extending 
the calculation of the cross section to higher missing energies yields the expected rise of 
high missing-momentum components in the cross section in comparison to the results near 
the Fermi energy. The contribution of various partial waves is studied demonstrating the 
increasing importance of higher /-values with increasing missing momentum. All these results 



are discussed in Sec. IV and a brief summary is presented in Sec. V. 



II. THE SINGLE-PARTICLE SPECTRAL FUNCTION 

The calculation of the cross section for exclusive (e, e'p) processes requires the knowledge 
of the hole spectral function which is defined in the following way 

>S(p,m„m^,p',ml,m^;£;) = | a^(p',m;,mO | *^-') (*^-' | a(p,m„mO | 

n 

5{E-{E^-Et')), (1) 

where the summation over n runs over the discrete excited states as well as over the con- 
tinuum of the (A-1) particle system, ^ is the ground state of the initial nucleus and 
a(p, ms,mT-) (a^(p', m^, m,-)) is the annihilation (creation) operator with the specified sp 
quantum numbers for momenta and third component of spin and isospin, respectively. The 
spectral function is diagonal in the third component of the isospin and ignoring the Coulomb 
interaction between the protons, the spectral functions for protons and neutrons are identi- 
cal for N=Z nuclei. Therefore in the following we have dropped the isospin quantum number 
m-j-. Note that the energy variable E in this definition of the spectral function refers to minus 
the excitation energy of state n in the A-1 particle system with respect to the ground-state 
energy {Eq ) of the nucleus with A nucleons. 

To proceed further in the calculations it is useful to introduce a partial wave decompo- 
sition which yields the spectral function for a nucleon in the sp basis with orbital angular 
momentum total angular momentum j, and momentum p 



SiAp,P; ^) = E (^0 I «k I ^n-') i^n-' I V. I ^0 ) S{E - {E^ - Ef)), (2) 

n 

where cipij{ap,ij) denotes the corresponding removal (addition) operator. The spectral func- 
tions for the various partial waves, Sij{p,p'; E), have been obtained from the imaginary part 
of the corresponding sp propagator gij{p,p';E). This Green's function solves the Dyson 
equation 
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gij{pi,P2] E) = glj\pi,p2; E) + J dps J dpiglj\pi,p3; E)A'Lij{ps,Pi; E)gij{p4,p2; E), (3) 

where g^^^ refers to a Hartree-Fock propagator and ASjj represents contributions to the 
real and imaginary parts of the irreducible self-energy, which go beyond the Hartree-Fock 
approximation of the nucleon self-energy used to derive g^^\ Although the evaluation of 
the self-energy as well as the solution of the Dyson equation has been discussed in detail in 
previous publications P,pi| we include here a brief summary of the relevant aspects of the 
method. 

A. Calculation of the nucleon self-energy 

The self-energy is evaluated in terms of a G-matrix which is obtained as a solution of 
the Bethe-Goldstone equation for nuclear matter choosing for the bare NN interaction the 
one-boson-exchange potential B defined by Machleidt (Ref. p3[, Table A. 2). The Bethe- 
Goldstone equation has been solved for a Fermi momentum kp = 1.4 fm~^ and starting 
energy —10 MeV. The choices for the density of nuclear matter and the starting energy 
are rather arbitrary. It turns out, however, that the calculation of the Hartree-Fock term 
(Fig. la) is not very sensitive to this choice Furthermore, we will correct this nuclear 
matter approximation by calculating the two-particle-one-hole (2plh) term displayed in 
Fig. lb directly for the finite system. This second-order correction, which assumes harmonic 
oscillator states for the occupied (hole) states and plane waves for the intermediate unbound 
particle states, incorporates the correct energy and density dependence characteristic of a 
finite nucleus G-matrix. To evaluate the diagrams in Fig. 1, we need matrix elements 
in a mixed representation of one particle in a bound harmonic oscillator while the other 

one can 



is in a plane wave state. Using vector bracket transformation coefficients 
transform matrix elements from the representation in coordinates of relative and center-of- 
mass momenta to the coordinates of sp momenta in the laboratory frame in which the two 
particle state is described by 
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\p1hj1P2hj2JT) (4) 

where pj, li and ji refer to momentum and angular momenta of particle i whereas J and 
T define the total angular momentum and isospin of the two-particle state. Performing an 
integration over one of the pi, one obtains a two-particle state in the mixed representation, 

/•OO 

\n1hj1p2l2j2JT) = dpiplRn^^ii{api) \p1l1j1p2l2j2JT) . (5) 
Jo 

Here Rni,ii stands for the radial oscillator function and the oscillator length a = 1.72 fm~^ 
has been chosen to have an appropriate description of the bound sp states in ^^O. Using the 
notation defined in Eqs. (^) and (^), our Hartree-Fock approximation for the self-energy is 
obtained in the momentum representation, 

^fuM^P'i) = 2(2-\l) ^ i2J+l){2T+l){p,hj,n2l2j2JT\G\p[hnn2l2j2JT). (6) 

y ' ' n2l2j2JT 

The summation over the oscillator quantum numbers is restricted to the states occupied in 
the independent particle model of ^^O. This Hartree-Fock part of the self-energy is real and 
does not depend on the energy. 

The terms of lowest order in G which give rise to an imaginary part in the self-energy 
are represented by the diagrams displayed in Figs. P) and |l]c, referring to intermediate 2plh 
and 2hlp states respectively. The 2plh contribution to the imaginary part is given by 

WSf'iPi,p[;E)= E E E jedkJK'dK{2J+l){2T + l) 

2{2ji + 1) n2l2j2 IL JJsST 

X {pihjin2l2j2JT\ G \klSJsKLT) {klSJsKLT\ G \p[hjin2l2j2JT) 

xrr6^E + e„,,,, -^-^-), (7) 

where the "experimental" sp energies en2i2j2 used for the hole states (—47 MeV, —21.8 



MeV, —15.7 MeV for s|, p| and p^ states, respectively), while the energies of the particle 
states are given in terms of the kinetic energy only. The plane waves associated with the 
particle states in the intermediate states are properly orthogonalized to the bound sp states 



following the techniques discussed by Borromeo et al. |2^. The 2hlp contribution to the 
imaginary part W^^^^{pi,p'i; E) can be calculated in a similar way (see also ||26|| ). 
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Our choice to assume pure kinetic energies for the particle states in calculating the 
imaginary parts of ly^P^'^ (Eq. (^) and VT^'^^p may not be very realistic for the excitation 
modes at low energy. Indeed a sizable imaginary part in ly^^^P is obtained only for energies 
E below —40 MeV. As we are mainly interested, however, in the effects of short-range 
correlations, which lead to excitations of particle states with high momentum, the choice 
seems to be appropriate. A different approach would be required to treat the coupling to 
the very low-lying 2plh and 2hlp states in an adequate way. Attempts at such a treatment 
can be found in Refs. |27|-^|6|. The 2plh contribution to the real part of the self-energy 
can be calculated from the imaginary part ly^P^"^ using a dispersion relation 



where V represents a principal value integral. A similar dispersion relation holds for y^'^^P 
and iy2hip_ 

Since the Hartree-Fock contribution has been calculated in terms of a nuclear matter 
G-matrix, it already contains 2plh terms of the kind displayed in Fig. |I]b. In order to avoid 
such an overcounting of the particle-particle ladder terms, we subtract from the real part of 
the self-energy a correction term (K), which just contains the 2plh contribution calculated 
in nuclear matter. Summing up the various contributions we obtain for the self-energy the 
following expressions 

S = S^^ + AS = + (y^plh _y^^ ^2hlp^ ^ |^p^2plh ^ p^2hlpj _ (g) 



B. Solution of the Dyson equation 

The next step is to solve the Dyson equation (|^) for the sp propagator. To this aim, we 
discretize the integrals in this equation by considering a complete basis within a spherical 
box of a radius -Rbox- The calculated observables are independent of the choice of -Rbox, if 
it is chosen to be around 15 fm or larger. A complete and orthonormal set of regular basis 
functions within this box is given by 
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^iijm{r) = {r\piljm) = Nu ji{pir) yijm{0, 0). (10) 

In this equation yijm represent the spherical harmonics including the spin degrees of freedom 
and ji denote the spherical Bessel functions for the discrete momenta pi which fulfill 

jiiPiRbo.) = 0. (11) 

Note that the basis functions defined for discrete values of the momentum pi within 
the box differ from the plane wave states defined in the continuum with the corresponding 
momentum just by the normalization constant, which is for the latter. This enables us 
to determine the matrix elements of the nucleon self-energy in the basis of Eq. (0) from 
the results presented in the preceding Subsection. 

As a first step we determine the Hartree-Fock approximation for the sp Green's function 
in the "box basis" . For that purpose the Hartree-Fock Hamiltonian is diagonalized 

max ^2 

E {P^\ ^^^n + \pn) {Pn\a),^ = ell • (12) 

Here and in the following the set of basis states in the box has been truncated by assuming 
an appropriate A^max- In the basis of Hartree-Fock states the Hartree-Fock propagator 
is diagonal and given by 

E - Cij ± 

where the sign in front of the infinitesimal imaginary quantity ir] is positive (negative) if e^^^ 
is above (below) the Fermi energy. With these ingredients one can solve the Dyson equation 
(0). One possibility is to determine first the socalled reducible self-energy, originating from 
an iteration of AS, by solving 

{a\ ^f{E) 1/3) = (a| AS,,(E) + ^ («| AS,,(i5) I7) gfir, E) (7I ^f{E) (14) 

7 

and obtain the propagator from 

gi,{a, E) = 6^,^ g\f{a- E) + gf^\a- E) {a\ J:]f{E) \(3) gl^\(3; E). (15) 
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Using this representation of the Green's function one can calculate the spectral function in 
the "box basis" from 

SfjiPm^Pn, ^) = ^ Im {Pm\a)ij 9ijia, E) {f3\pn)i^ ■ (16) 

For energies E below the lowest sp energy of a given Hartree-Fock state (with Ij) this spec- 
tral function is different from zero only due to the imaginary part in Y7°^. This contribution 
involves the coupling to the continuum of 2hlp states and is therefore nonvanishing only 
for energies at which the corresponding irreducible self-energy AS has a non-zero imaginary 
part. Besides this continuum contribution, the hole spectral function also receives contri- 
butions from the quasihole states 0. The energies and wave functions of these quasihole 
states can be determined by diagonalizing the Hartree-Fock Hamiltonian plus AS in the 
"box basis" 

ax 

E (P^l ^^^n + + AS,,(E = e%) b„) (pJT),^. = e% (p,|T),^. . (17) 

n=l 

Since in the present work AS only contains a sizable imaginary part for energies E below e^^, 
the energies of the quasihole states are real and the continuum contribution to the spectral 
function is separated in energy from the quasihole contribution. The quasihole contribution 
to the hole spectral function is given by 

SfijiPm^Pn, E) = ZrijiPmWlj {^\Pn)ijS{E - efij), (18) 

with the spectroscopic factor for the quasihole state given by [^] 



/ g(T|AS,,(E)|T) 



(19) 



Finally, the continuum contribution of Eq. (0) and the quasihole parts of Eq. (|T8D, which 
are obtained in the basis of box states, can be added and renormalized to obtain the spectral 
function in the continuum representation at the momenta defined by Eq. (Ill]) 



2 1- 

Slj{p^,Pn]E) = - ^ {Sij {p,n,Pn;E) Sflj {Pm,Pn]E)). (20) 
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It is useful to have a separable representation of the spectral function in momentum space. 
For a given energy, the spectral function in the box is represented by a matrix in momentum 
space; after diagonalizing this matrix one obtains 

^ max 

Sij{pm,Pn;E) = J2 (piiPm) MPn) (21) 

i 

where Sij{i) are the eigenvalues and 0, are the corresponding eigenf unctions. In all cases 
considered here, it is enough to consider the first 5 or 6 largest eigenvalues in Eq. (^l]) for 
an accurate representation of the spectral function. These eigenfunctions are in principle 
sp overlap functions (see discussion after Eq. (129|) below). They can be thought of as the 
natural orbits at a given energy. In fact, if the diagonalization is performed after integrating 
over the energy E one would precisely obtain the natural orbits associated with the one-body 
density matrix and the eigenvalues Sij{i) would be the natural occupation numbers |[T0|| . 



III. GENERAL FORMALISM OF DWIA 

For the scattering of an ultrarelativistic electron with initial (final) momentum pe (Pg), 
while a nucleon is ejected with final momentum p^, the differential cross section in the 



one-photon exchange approximation reads |119| , |2T 

da 1 



S Lx,yWx,y, (22) 



where = — and q = Po — Pc? ^ = Pc~P'e momentum and energy transferred 

to the target nucleus, respectively. The quantities Lx^\',W\^x' (usually referred to as the 
lepton and hadron tensors, respectively) are expressed in the basis of unit vectors 

eo = (l,0,0,0), 
e±i = (0,TV^,-vTi,0), (23) 

which define the longitudinal (0) and transverse (±1) components of the nuclear response 
with respect to the polarization of the exchanged virtual photon. The components of the 
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lepton tensor depend only on the electron kinematics, while Wx^y depend on q, Lj,p'j^, cos 7 
Pn ■ q/Pn'?5 ^iid the angle a between the (p^, q) plane and the electron scattering plane. 



The hadron tensor is defined as |19,|2T|,|33| 



w,,y = E f ^A(q)^v(q) s (^i - ^f) , (24) 

i Ji 

i.e. it involves the average over initial states and the sum over the final undetected states 
(compatible with energy-momentum conservation) of bilinear products of the scattering 
amplitude JA(q)- 

This basic ingredient of the calculation is defined as 

J,(q)= Jdre^^-'^{^t\J,-e',K), (25) 

where the matrix element of the nuclear charge-current density operator is taken between 
the initial, I^I/q ), and the final, ), nuclear states. A natural choice for is suggested by 
the experimental conditions of the reaction selecting a final state, which behaves asymptot- 
ically as a knocked out nucleon with momentum p'j^ and a residual nucleus in a well-defined 
state I^^^H^)) with energy E and quantum numbers n. By projecting this specific channel 
out of the entire Hilbert space, it is possible to rewrite Eq. (^) in a one-body representation 
(in momentum space and omitting spin degrees of freedom for simplicity) as [l^ 



J,(q) = /dp xj,7^;(p + q) Jfip, q) " </'^„(p)[5n(i?)]^ (26) 

provided that is substituted by an appropriate effective one-body charge-current density 
operator J^^, which guarantees the orthogonality between and ) besides taking into 
account effects due to truncation of the Hilbert space. Actually, the orthogonality defect 
is negligible in the standard kinematics for (e, e'p) reactions and in DWIA J^^ is usually 
replaced by a simple one-body current operator [p!8| , pO| , pT . 
The functions 



4v^LiP) = {^n-\EMp)K) (27) 
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describe the overlap between the residual state ^(i?)) and the hole produced in |\I/q ) 
and I'^f), respectively, by removing a particle with momentum p. Both 4'EnyXp'^En 
eigenfunctions of a Feshbach-like nonlocal energy- dependent Hamiltonian referred to the 



residual nucleus, belonging to the eigenvalues E and E + uj, respectively |T^. The norm of 
is 1 and Sn{E) is the spectroscopic factor associated with the removal process, i.e. it is 
the probability that the residual nucleus can indeed be conceived as a hole produced in the 
target nucleus. The dependence of Xp' ^En ^'n hidden in the asymptotic state ) and 
the boundary conditions are those of an incoming wave. 

Because of the complexity of the eigenvalue problem in the continuum, the Feshbach 
hamiltonian is usually replaced by a phenomenological local optical potential V{r) of the 
Woods-Saxon form with complex central and spin-orbit components. It simulates the mean- 
field interaction between the residual nucleus and the emitted nucleon with energy- dependent 
parameters determined through a best fit of elastic nucleon-nucleus scattering data including 
cross section and polarizations. Then, X^/^En ~ ^p~^ expanded in partial waves and 
a Schrodinger equation including V^(r) is solved for each component up to a maximum 



angular momentum satisfying a p'^-dependent convergency criterion . The nonlocality 
of the original Feshbach hamiltonian is taken into account by multiplying the optical-model 
solution by the appropriate Percy factor |^ . 



After summing over the undetected final states with quantum numbers n of the residual 
nucleus, the hadron tensor Wx^y in momentum space becomes 

W^A,A' ~ E / dpdp' xil^*ip + cOJ.iP^ q) ■ 0^„(p)0k(p')^n(i?) 

n 

^ J dpdp' xil^*{p + q)i^(p, q) ■ ^(p, P; E) 

^](p',q)-e^;x^7^(p' + q), (28) 



where 



5(p, p'; E)=Y: 5„(E)0^„(p')0^„(p) (29) 

n 
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is the hole spectral function defined in Eq. (|l]). Notice that the spin and isospin indices 
have been omitted for simplicity and the summation over n is over the different partial wave 
contributions which are present at a given energy E. This sum should not be confused with 
the separable representation (Eq. (|2l|)) of the partial wave contributions to the spectral 
function Sij{p,p', E) defined in Eq. (^). Each /j-contribution, coming from either quasi-hole 
states (if E is the correct excitation energy) or from states which are usually unoccupied in 
the standard shell model, can be separately computed, so that the total hadron tensor will 
look like 

W^A,A' = E <A' ■ (30) 

Experimental data for the (e, e'p) reaction are usually collected as ratios between the 
measured cross section and Ka^^, where is a suitable kinematical factor and cTeN is 
the elementary (half off-shell) electron-nucleon cross section. In this way the information 
contained in the five-fold differential cross section is reduced to a two-fold function of the 
missing energy E^ = uo — Tp'^ — E^ (Tpj^ is the kinetic energy of the emitted nucleon and E^ is 
the excitation energy of the residual nucleus) and of the missing momentum Pm = — q . 
Therefore, in the following Section results will be presented under the form of the socalled 



reduced cross section [21 



, , dfj 1 . 
dp^^dpj^ iTcTeN 



IV. RESULTS 

In this Section we will discuss results for the reduced cross section defined in Eq. ( |3TD 
for (e, e'p) reactions on ^^O leading both to discrete bound states of the residual nucleus ^^N 
and to states in the continuum at higher missing energy. Distortion of electron and proton 



waves has been taken into account through the effective momentum approximation |35 



and through the optical potential derived from the Schwandt parametrization [BB] (see Tab. 
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Ill in Ref. 0), respectively. All results presented here have been obtained using the CCl 
prescription for the half off-shell elementary electron-proton scattering amplitude in 
analogy with what has been commonly done in the analysis of the experimental data. We 
also employed the nonrelativistic description for this amplitude |^ to be consistent with the 
nonrelativistic calculation of the five-fold differential cross section. In parallel kinematics, 
where most of the experimental data are available, this choice does not produce very different 
results with respect to the former, and, therefore, will not be considered in the following. 

A. Quasihole states 

In Fig. H the experimental results for the transition to the ground state of ^^N are dis- 
played as a function of the missing momentum p^. These data points have been collected at 
NIKHEF choosing the socalled parallel kinematics [Q, where the direction of the momentum 
of the outgoing proton, p^, has been fixed to be parallel to the momentum transfer q. In 
order to minimize the effects of the energy dependence of the optical potential describing 
the FSI, the data points have been collected at a constant kinetic energy of 90 MeV in the 
center-of-mass system of the emitted proton and the residual nucleus. Consequently, since 
the momentum of the ejected particle is also fixed and 

Pm = IpnI - kl> (32) 

the missing momentum can be modified by collecting data at various momenta q transferred 
from the scattered electron. 

The experimental data points for this reduced cross section are compared to the pre- 
dictions of the calculations discussed above. The quasihole part of the spectral function 
for the p| partial wave represents the relevant piece of the nuclear structure calculation for 
the proton knockout reaction leading to the ground state of ^^N. Using the quasihole part 
of the spectral function as discussed above (see Eq. (p!8|)) but adjusting the spectroscopic 
factor for the quasihole state contribution Zq i to fit the experimental data, we obtain the 
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solid line of Fig. |^. Comparing this result with the experimental data one finds that the 
calculated spectral function reproduces the shape of the reduced cross section as a function 
of the missing momentum very well. The absolute value for the reduced cross section can 
only be reproduced by assuming a spectroscopic factor Zq^i = 0.644, a value considerably 
below the one of 0.89 calculated from Eq. (|I9|) 0]. The phenomenological Woods-Saxon 
wave functions adjusted to fit the shape of the reduced cross section require spectroscopic 
factors ranging from 0.61 to 0.64 for the lowest Op| state and from 0.50 to 0.59 for the 
Op| state, respectively, depending upon the choice of the optical potential for the outgoing 
proton [^. The fact that the calculated spectroscopic factor is larger than the one adjusted 
to the experimental data may be explained by the observation that the calculations only 
refiects the depletion of the quasihole occupation due to short-range correlations. Further 
depletion and fragmentation should arise from long-range correlations due to collective exci- 



tations at low energies ||6|,|3T||. Other explanations for this discrepancy could be the need for 
improving the description of spurious center-of-mass motion [^,0] or a different treatment 
of FSI in terms of a relativistic model for the optical potential ||40|| . 



In order to visualize the effects of FSI, Fig. |^ also displays the results obtained for the 
quasihole contribution to the spectral function (with the same spectroscopic factor Z^^i = 
0.644 as before, for sake of consistency) but ignoring the effects of the optical potential. 
In this socalled Plane- Wave Impulse Approximation (PWIA) the reduced cross section as a 
function of the missing momentum is identical to the spectral function at the missing energy 
of the considered Op| state, or, better, to the momentum distribution of the peak observed 
at this missing energy with the quantum numbers of the ground state of ^^N. Therefore, 
the difference between the solid and the dashed line in Fig. ^ corresponds to the difference 
between the reduced cross section defined in Eq. ( pi]) and the momentum distribution for 
the ground state of ^^N. In other words, it illustrates the effect of all the ingredients entering 
the present theoretical description of the (e, e'p) reaction, which are not contained in the 
calculation of the spectral function. In particular, the real part of the optical potential 
yields a reduction of the momentum of the outgoing proton p'^. According to Eq. (|32D , this 
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implies in parallel kinematics a redistribution of the strength towards smaller values of the 
missing momentum and makes it possible to reproduce the observed asymmetry of the data 
around Pm = 0. This feature cannot be obtained in PWIA (dashed line), where the results 
are symmetric around Pm = due to the cylindrical symmetry of the hadron tensor y 
around the direction of q when FSI are switched off (for a general review see Ref. and 



references therein). The imaginary part of the optical potential describes the absorption 
of the proton flux due to coherent inelastic rescatterings, which produces the well known 
quenching with respect to the PWIA result. 

As a second example for the reduced cross section in (e, e'p) reactions on ^^O leading to 
bound states of the residual nucleus, we present in Fig. |^ the data for the | state of ^^N at 
an excitation energy of —6.32 MeV. Also in this case the experimental data are reproduced 
very well if we adjust the spectroscopic factor for the corresponding quasihole part in the 
spectral function to Zq^s = 0.537. The discrepancy with the calculated spectroscopic factor 
(0.914) is even larger for this partial wave than it is for the p| state. A large part of this 
discrepancy can be attributed to the long-range correlations, which are not accounted for in 
the present study. Note, that in the experimental data three | states are observed in ^^N 
at low excitation energies. Long-range correlations yield a splitting such that 86% of the 
total strength going to these three states is contained in the experimental data displayed 
in Fig. ^ This splitting is not observed in the theoretical calculations. If one divides the 
adjusted spectroscopic factor Z^^a by 0.86 to account for the splitting of the experimental 
strength, one obtains a value of 0.624 which is close to the total spectroscopic factor adjusted 
to describe the knockout of a proton from state. 

Figure ^ also contains the results for the reduced cross section derived by substituting 
the overlap [5'„(-E')]5 0^^ in Eq. ( p6D with the variational wave function of Pieper et al. 0, 
who employed the Argonne potential for the NN interaction |^ . Also in this case the shape 
of the experimental data is globally reproduced with a slightly better agreement for small 
negative values of pm but with a clear underestimation at larger p^. The overall quality of 
the fit is somewhat worse than for the Green's function approach and the required adjusted 



spectroscopic factor is Zq^^ = 0.459, even below the value of 0.537 needed in the present 



2 

calculation. 



The analysis of the reduced cross section has been extended to higher missing momenta 



by experiments performed at the MAMI accelerator in Mainz adopting different kine- 
matical conditions than the parallel kinematics. Using the same spectroscopic factors for the 
p| and the p| partial waves, which were adjusted to the NIKHEF data above, the results 
of our calculations agree quite well also with these MAMI data, as displayed in Fig. ^. Al- 
though the calculation is somewhat below the data at high missing momentum, one should 
keep in mind that the corresponding difference in sp strength is only an extremely tiny 
fraction of the 10% of the protons which are expected to be associated with high momenta 
due to short-range correlations [BHTO . 



B. The contribution of the continuum 

/,From theoretical studies it is known that an enhancement of the high-momentum com- 
ponents due to short-range NN correlations does not show up in knockout experiments 
leading to states of low excitation energy in the (A-1) nucleus, but should be seen at higher 
missing energies, which correspond to large excitation energies in the residual nucleus. A 
careful analysis of such reactions leading to final states above the threshold for two-nucleon 
emission, however, is much more involved. For example, a description of the electromag- 
netic vertex beyond the impulse approximation is needed and two-body current operators 
must be adopted which are consistent with the contributions included in the spectral func- 
tion. Moreover, the possible further fragmentation of the (A-1) residual system requires, 
for a realistic description of FSI, a coupled-channel formalism with many open channels. 
Calculations based on the optical potential are not satisfactory at such missing energies, 
because inelastic rescatterings and multi-step processes will add and remove strength from 
this particular channel. 

Nevertheless, it should be of interest to analyze the predictions of the present approach 



19 



at such missing energies. First of all, because it represents the first realistic attempt of 
a complete calculation of the single-particle channel leading to the final proton emission, 
including intermediate states above the Fermi level up to / = 4; therefore, it represents a 
realistic estimate of the relative size of this specific channel. Secondly, because information 
on the shape of the reduced cross section as a function of the missing momentum or on the 
relative contribution of various partial waves could yield reliable results even at these missing 
energies. Due to the problems mentioned above, no reliable description of the absolute value 
of the reduced cross section can be reached in this framework. 

In order to demonstrate the energy dependence of the spectral function and its effect on 
the cross section, we have calculated the reduced cross section for the excitation of | states 
at = —63 MeV. For these studies we considered the socalled perpendicular kinematics, 
where the energy of the emitted proton is kept fixed at 90 MeV as well as the momentum 
transfer at g ~ 420 MeV/c (equal to the outgoing proton momentum). The same opti- 
cal potential as in Figs. 0, ^ can be adopted to describe FSI and the missing momentum 
distribution is obtained by varying the angle between and q. For a spectral function 
normalized to unity (as the absolute result for the cross section is not reliable), the reduced 
cross section is represented by the solid line in Fig. ^ If, however, we replace the spectral 
function derived from the continuum contribution in Eq. (^OD by the one derived for the 
I quasihole state at its proper missing energy (but now in the same kind of perpendicular 
kinematics and normalized to 1) we obtain the dashed line. A comparison of these two 
calculations demonstrates the enhancement of the high-momentum components in the spec- 
tral function leading to final states at large excitation energies. Note that the cross section 
derived from the appropriate spectral function is about two orders of magnitude larger at 
Pm ~ 500 MeV/c than the one derived from the spectral function at the quasihole energy. 

The discussion so far is of course somewhat academic since it will be difficult to perform 
a decomposition of the continuum contribution to the reduced cross section in terms of the 
quantum numbers for angular momentum and parity of the state for the residual system. 
Therefore we display in Figs. ^ and the contributions to the total reduced cross section of 
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the various partial waves associated to states above the Fermi level and usually unoccupied in 
the standard shell model. From Fig. ^we can furthermore see that the relative importance of 
the various partial waves changes with the missing momentum, emphasizing the contribution 
of higher angular momenta at increasing p^. This feature can be observed even better in 
Fig. 1^, where the percentage of each relative contribution to the total reduced cross section 
is displayed as a function of the missing momentum. For each orbital angular momentum 
we obtain a "window" in where its contribution shows a maximum as compared to other 
partial waves. 



V. CONCLUSIONS 

In the present paper the consequences of the presence of high-momentum components in 
the ^^O ground state have been explored in the calculation of the (e, e'p) cross section within 
the formalism for the DWIA developed in Refs. p!7|-p2|. The spectral functions have been 
calculated for the ^^O system itself, by employing the techniques developed and discussed 
26|,P-p^. At low missing energies, the description of the missing momentum dependence 



m 



of the p| and p| quasihole states compares favorably with the experimental data obtained 
at NIKHEF P] and at the MAMI facility in Mainz |15|. The difference between theory and 



experiment at high missing momenta can at most account for a very tiny fraction of the 
sp strength which is predicted to be present at these momenta [p|-|lO[]. A comparison with 
the PWIA result clarifies the influence of FSI in parallel kinematics. We also compare our 
results for the p| quasihole state with the results obtained in Ref. for the Argonne NN 
interaction. While the shape of the cross sections is nicely described by our results, the 
associated spectroscopic factors are overestimated substantially. Although a large fraction 
of this discrepancy can be ascribed to the influence of long-range correlations PJ31[] , which 



are outside the scope of the present work, a discrepancy may still remain although it has 
been suggested that a correct treatment of the center-of-mass motion |^ may fill this gap. 
As discussed previously for nuclear matter (see e.g. and emphasized in P|-pi]|] for 
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finite nuclei, the admixture of high- momentum components in the nuclear ground state can 
only be explored by considering high missing energies in the (e, e'p) reaction. Although 
other processes may contribute to the cross section at these energies, we have demonstrated 
in this paper that the expected emergence of high missing momentum components in the 
cross section is indeed obtained and yields substantially larger cross sections than the cor- 
responding outcome for the quasihole states. As a result, we conclude that the presence of 
high-momentum components leads to a detectable cross section at high missing energy. In 
addition, we observe that it is important to include orbital angular momenta at least up 
to / = 4 in the spectral function in order to account for all the high missing momentum 
components up to about 600 MeV/c. A clear window for the dominant contribution of each 
/-value as a function of missing momentum is also established. This feature may help to 
analyze experimental data at these high missing energies. 
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FIGURES 

FIG. 1. Graphical representation of the Hartree-Fock (a), the two-particle-one-hole (2plh,b) 
and the two-hole-one-particle contribution (2hlp,c) to the self-energy of the nucleon. 

FIG. 2. Reduced cross section for the ^^0(e, e'p)^^Ngs reaction in parallel kinematics. Results 
with (solid line) and without (dashed line) inclusion of the FSI are compared to the experimental 
data [|[. A spectroscopic factor of 0.644 has been employed in displaying the results for the 
calculations involving the spectral function. 

FIG. 3. Reduced cross section for the ^^0{e, e'p) reaction in parallel kinematics leading to the 
I state at —6.32 MeV of the residual nucleus ^^N. Results of the present Green's function approach 
(solid line) are compared to those obtained in the variational calculation of 0] (dashed line) and 
the experimental data |p. A spectroscopic factor of 0.537 was required for the Green's function 

approach, while Z 3 = 0.459 has been used to adjust the results of the variational calculation. 

P2 

FIG. 4. Reduced cross section for the ^^0{e, e'p) reaction leading to the ground and the | 
states of ^^N in the kinematical conditions considered in the experiment of The calculations 
were performed using the same spectral functions as discussed for Figs. 2 and 3. 

FIG. 5. Reduced cross section for the ^^0{e,e'p) reaction in perpendicular kinematics for the 
excitation of | states at Ey^ = —63 MeV (solid line) and —6.32 MeV (dashed line). 

FIG. 6. Contributions of various partial waves to the reduced cross section for the ^^0{e, e'p) 
reaction in the same conditions as for the solid line in Fig. 5. 

FIG. 7. Relative importance of various partial waves to the reduced cross section for the 
^^0{e, e'p) reaction in the same conditions as in Fig. 6. 
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